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Abstract. It is the aim of the paper to present a new point of view on ro- 
tational elasticity in a nonlinear setting using orthogonal matrices. The pro- 
posed model, in the linear approximation, reduces to the well known equi- 
librium equations of static linear elasticity. An appropriate kinetic energy is 
identified and we present a dynamical model of rotational elasticity. The prop- 
agation of elastic waves in such a medium is studied and we find two classes of 
waves, transversal rotational waves and longitudinal rotational waves, both of 
which are solutions of the nonlinear partial differential equations. For certain 
parameter choices, the transversal wave velocity can be greater than the lon- 
gitudinal wave velocity. Moreover, parameter ranges are identified where the 
model describes an auxetic material. However, in all cases the potentional en- 
ergy functional is positive definite. We present one possible way of visualising 
transversal rotational waves. 



1. Introduction 

We are developing a theory of rotational elasticity based on using orthogonal 
matrices as the basic dynamical variables. Let us consider a three-dimensional 
clastic continuum, occupying an open connected set f2 C R 3 or the whole of K 3 , 
whose material points can only experience rotations and no displacements. Our 
Euclidean space is assumed to be equipped with Cartesian coordinates x = (x, y, z). 
To every point x G R 3 we attach an orthonormal basis {ei, 62,63}. In the initial 
state these basis vectors are assumed to be aligned with our Cartesian coordinate 
axes, i.e. (ej)j = 5ij, but will become functions of x and time t once the continuum 
is deformed. 

Models of this type have in fact a long tradition. One such model has first been 
introduced by MacCullagh in 1839, see [27]. He noted that is was not possible to 
describe optical phenomena by comparing the aether with an ordinary elastic solid. 
He thus introduced a new type of medium whose potential energy depended only on 
rotations. Similar models have been investigated by the Cosserat brothers which 
resulted in an extended framework of elasticity [5], often referred to as Cosserat 
elasticity. The main difference between Cosserat elasticity and classical elasticity 
is the assumed independence of displacements and rotations, often referred to as 
microrotations. When formulating an energy functional of classical elasticity, one 
assumes the integrand not to depend on the derivatives of rotation. It is well 
known that rotations about different axes do not commute, and thus we expect an 
inherently nonlinear theory. Variants of the theory of Cosserat elasticity appear 
under various names in modern applied mathematics literature such as oriented 
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medium, asymmetric elasticity, micropolar elasticity to mention a few popular ones, 
see e.g. [HI ESI [El El [231 18] . 

This field is also related to the theory of granular media, ferromagnetic materials, 
cracked media and liquid crystals to mention a few of them. There has also been an 
interest in the theory of Cosserat elasticity from a more theoretical physics point 
of view, see [241 111 Ull Hi]- 

In this work we are only considering rotational deformations of the continuum 
and neglect displacements. We assume that every material point can experiences 
an independent rotation e* n- 0(x)e l , i = 1,2,3, where O(o;) denotes the field 
of orthogonal matrices, O t O = OO t = I where O t denotes the transpose and I 
the identity matrix. We consider orthogonal matrices with detO = +1. The most 
general potential energy functional of such a static continuum is thus taken to be 



(1.1) V := [ f(x,0,dO,d 2 0,...)dxdydz 

Jn 



where dO (d 2 0, . . .) denotes the first (second, . . .) partial derivatives of O. The 
matrix- function 0(a;) is the unknown quantity, our dynamical variable. 

2. Statement of the static problem 

2.1. Basic equations of the static model. The general framework (|1.1|) of our 

rotational elasticity model is very wide. Hence, as a first step we identify a suitable 
subclass of variational functionals which is of mechanical and physical interest, and 
which is invariant under rigid rotations. This means the mapping O M> OO, where 
O is a constant orthogonal matrix, leaves the variational functional invariant. 
Possible deformations of our elastic medium are characterised by the quantity 

where we sum over repeated indices. In the context of Cosserat elasticity, this 
quantity is sometimes referred to as contortion jTH [17] ■ Note that K, 3 t is, by 
definition, skew-symmetric in the first and third index, K^-fc = — K^. This follows 
from the fact that 9(OO t ) = 0. Moreover, the quantity K is invariant under rigid 
rotations. Since we are working in R 3 it seems natural to contract K, 3J t with the 
Levi-Civita symbol to construct a 3 x 3 matrix to analyse the deformations, we 
define 

Since the matrix A has no a priori symmetries, it is useful to decompose it into its 
three irreducible components, a trace part, a skew-symmetric part and a trace-free 
symmetric piece, respectively defined by 



mn 


— -A 1 


trace 


A (2) 

mn 




skew — symmetric 


A (3) 

mn 


A _ A(l) _ A(2) 

■— -n-mn ^-mn ^mn 


trace — free symmetric 



(3) 

Note that the definition for A mn can also be written in the form 



A^ 3 ) — -(A +A "I--A-T 
which makes explicit that this piece is indeed trace-free and symmetric. 
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This trace-free symmetric piece is closely related to the so-called Q-tensor of 
the Landau-de Gennes theory of uniaxial nematic liquid crystals, see [2j[2Q]. The 
Q-tensor is symmetric, traceless and in addition has two equal eigenvalues. The 
elastic energy density in this theory is based on elastic invariants constructed from 
the irreducible components of the first partial derivatives of Q. In our model this 
corresponds to d 2 0, the second partial derivatives of O. Therefore, our general 
potential energy is sufficiently wide to cover a large variety of materials. 

A natural starting point would be to consider a functional of the type 

V x = [ [dllA^^f-l-callA^f + callA^II 2 ! dxdydz 
Jn L J 

where ci, c 2 and C3 are positive constants which are naturally called the elastic 
moduli. A similar functional was considered in [24j to formulate a unified field 
theory for mesons. Moreover, we denote the norm by ||AW|| 2 = Am„Am n where 
i = 1,2,3 and no summation over i. Note that a similar functional has been 
considered in the context of Cosserat elasticity [T3]. The three elastic moduli are 
not fully independent when one works in the whole of R 3 because of the following 
identity. One can show that 

(2.1) -2||A«|| 2 - ||A( 2 )|| 2 + IIA^f = AeijkdiAfl 

This identity has its roots in differential geometry. Of course, when f2 has a bound- 
ary, the above argument fails and the identity (12.11) leads to the appearance of 
boundary terms. One can regard K as an affine connection on a globally flat Rie- 
mannian manifold with curvature and torsiorQ. The definition of the Riemann 
curvature tensor involves derivatives of the connection and contraction of the con- 
nection with itself. The vanishing of the Riemann curvature tensor is then equiva- 
lent to (|2.1|) . Note that the vanishing of the Riemann curvature tensor is also often 
used in classical elasticity since one assumes that the distorted medium remains 
flat, thus implying certain compatibility conditions. 

Next we can rewrite || A^ 3 ' || 2 in the energy functional V\ in terms of the other two 
terms and a surface term which will not alter the resulting equations of motions. 
This results in the functional 



(2.2) V 2 



Cl ||A( 1 )|| 2 + c 2 ||A( 2 )|| 2 



dx dy dz 



where C\ = C\ + 2c$ and c 2 = c 2 + c 3 . 



2.2. Basic equations of the linearised model. It is known that every orthog- 
onal matrix O can be written as the exponential of a skew-symmetric matrix. It 
turns out to be convenient to analyse the nonlinear equations and their relation with 
linear elasticity by writing O = exp(it*). We denote by it* the a skew-symmetric 



1 Alternatively, one can work with the torsion tensor T directly [3], which is defined by T^j.; := 
Kjkl — Kjlk- This tensor can also be contracted with the Levi-Civita symbol to construct a 

3x3 matrix to be decomposed into irreducible pieces. Both formulations are equivalent and can 

* * * 

easily be translated into each other by noting that AM = -T ax , A< 2 ) = -2T voc , A< 3 ) = 2T ten , 
||6|| = V2\\oj\\ and identifying Cl = c ax , 4c 2 = c voc , 4c 3 = c tcn , 2p = c kin . 
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matrix dual to the vector it. This means (u*)ik — £ijkUj or, written out explicitly, 

— u z Uy 
u z — u a 

~Uy U X 

In other words, rather than using elements of the group O 6 SO (3) as the dynamical 
or state variables, we are using elements of its Lie algebra u* £ so (3), see [14] where 
this approach is used in linear Cosserat elasticity. 

Thus we can expand the orthogonal matrix in u and obtain 

(2.3) = exp(«*)=I + «* + 0(||M|| 2 ) 

where ||«|| 2 denotes the standard vector norm ||m|| 2 = u\ + u 2 + u 2 . It should be 
noted that there exits a simple explicit representation for the matrix O which is 
sometimes referred to as Rodrigues' formula 

exp(u*) = I + sin(||n||)— + (1 - C os(||u||))^-. 

HI HI 

Let us now assume u to be small and let us expand the terms in the func- 
tional (|2.2p up to terms 0(||u|| 3 ). This gives 

||AW|| a = |||divu|| 2 + 0(H| 3 ) 
||A( 2 )|| 2 =2||curlM|| 2 + 0(||w|| 3 ). 
Therefore, in the linear approximation the energy functional (12.2[) becomes 

(2.4) V 3 = I a\\divu\\ 2 + f3\\cur\u\\ 2 dxdydz 

where we denoted a — 4ci/3 and j3 = 2c2- 

We can also expand the term ||A^ 3 ^|| 2 in terms of the vector u. However, it is 
more useful to consider the right-hand side of (|2.ip which is given by 

(2.5) SijkdiA^ = d m u n d n u m - ||divM|| 2 + 0(||m|| 3 ). 

Variation of the basic linear functional (I2.4[) with respect to our independent 
(dynamical) variable, the components of u* , yields 

f ( dV 3 dV 3 \ 

SV3 = L ld< b SKb + w^) 5{dc<b) ) dxdydz ' 

which after integration by parts and using Gauss' theorem gives 

6V * = I {t^T ~ 9 ° am^* ^ ) S <b dx dy dz + f n ^ S(d c u* ab )dS. 
Jn\du* ab d(d c u* ab )J J dn d(d c u* ab ) 

To simplify matters, we assume Dirichlet boundary conditions and henceforth ne- 
glect the boundary term. The gives us the Euler-Lagrange equations 

a grad(div u) — ft curl curl u = 0. 

This system of equations has already been known by Cauchy and led to modifica- 
tions due to MacCullagh and also Neumann, we refer the reader to |27j . 

The vector identity curl curl u = grad(divtt) — Am enables us to rewrite the 
Euler-Lagrange equations in the form 

(2.6) ft Au + (a - ft) grad(div it) = 
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Comparing with the well-known equilibrium equations of linear elasticity /i Am + 
(/i + A) grad(div u) — 0, we have to choose ft = H and a = A + 2fi to get a corre- 
sponding theory. Therefore, we have established an analogy between the linearised 
theory of the nonlinear rotational elasticity model and classical linear elasticity. 
Therefore, one can now argue that the energy functional (|2.2p is in fact a 'natural' 
nonlinear generalisation of the functional of linear elasticity. Note that the equilib- 
rium equations of linear elasticity are usually derived from a different functional, 
see [HI [221 [25]. If we denote = (diUj + djUi)/2, then the energy functional is 
given by 



It should be noted here that parts of the second term of the functional can be found 
in the right-hand side of the derivative term (12.51) . Another interesting fact of linear 
elasticity is that not all the irreducible pieces of the matrix diUj are regarded as 
the fundamental building blocks of the potential energy functional. The skew- 
symmetric part, namely the rotations, are in general not part of the potential 
energy. It is the aforementioned Cosserat theory of elasticity which incorporates 
local rotations of points. 

The functional (|2.7[) is obviously positive definite if the elastic moduli A and fi 
are assumed to be positive. Similarly, the functional (|2.4p is also positive definite 
provided that a and (3 are positive. Therefore, when comparing the nonlinear 
functional (|2.2|) . which is positive definite provided c$ > for i = 1, 2, 3, with the 
linear elasticity functional (|2.7|) . we are led to identify 



This immediately suggests that Poisson's ratio for our rotational medium may 
not be between —1 and 1/2 and, moreover, that the transverse wave velocity may 
be greater than the longitudinal wave velocity. By straightforward comparison we 
find for Poisson's ratio a and Young's modulus E, respectively 



Thus, by choosing the constants c; > appropriately, it is possible to define ma- 
terials with negative Poisson's ratio, also known as auxetic materials, see [16]. For 
ordinary materials with positive Young's modulus E > and positive Poisson's 
ratio in the range a 6 (0, 1/2), our positive constants have to satisfy c\ > 3c2 + C3. 
An auxetic material with a £ (—1, 0) is realised when 2c2 < c\ < 3c2 + C3. 

Moreover, the constants Ci > can also be chosen so that the analogous linearised 
functional (|2.7p is not positive definite while the nonlinear theory in fact is. This 
fact is important when stability is investigated. 

2.3. Coercivity and convexity. The equilibrium equations of our rotational elas- 
ticity theory are derived from an energy functional. Therefore, the nonlinear model 
is automatically a quasilinear, second order PDE system in divergence form [12) . 
A natural topic to address is the convexity of the potential energy [1]. For the 
linearised model (|2.4[) one can immediately verify that the functional is coercive, 



(2.7) 






A ci — 3c 2 — c 3 



E 



2(A + fi) 2ci - 3c 2 + c 3 ' 
/i(3A + 2/i) _ 12(ci-2c 2 )(c 2 + c 3 ) 
A + 2ci - 3c 2 + c 3 
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independent of u and convex in du. Thus, existence and uniqueness of the lin- 
earised Dirichlet boundary value problem are established in appropriate function 
spaces. 

Since our basic variables of the nonlinear model are matrices, convexity of the 
functional will also depend on the convexity of the matrix exponential. The defi- 
nition of a matrix valued convex function was introduced in |15j and is analogous 
to that of a scalar function. It follows that the matrix exponential is not convex 
in general. Our starting point functional is constructed to contain only sums of 
squares, so it is clearly bounded from below by zero. 

However, coercivity of the functional remains a difficult task to be addressed 
since it does not seem to be true in general. By definition we have || exp(u*)|| 2 = 3, 
however, the norm ||u|| is apriori unconstrained. 

2.4. Identifying kinetic energy. So far our model of rotational elasticity has 
been static. The next logical step is to introduce kinetic energy into the energy 
functional (|2.2p . To do this, we add a term of the form ||0|| 2 to the energy func- 
tional. To be consistent with our previous notation, we have 

(2.8) ||6|| 2 = trace(66 T ) 

which is invariant under rigid rotations. If we now linearise this form of kinetic 
energy according to Eq. (|2.3[) . we find 

(2.9) ||6|| 2 = 2\\u\\ 2 + 0{u 3 ). 

Note that one regards ||0|| as angular speed (modulus of the angular velocity vector) 
and /?||0|| 2 as rotational energy or angular kinetic energy, where p is the energy 
density of the medium. 

2.5. Statement of the full problem. After identifying the kinetic energy, we 
can now formulate the complete variational functional of rotational elasticity in 
our setting. 

Let us find a skew-symmetric matrix u* such that O = exp(u*) extremises the 
variational functional or action 



(2.10) V= [ kllAWf + callAWf + callAWf-pHbl 



dx dy dz dt 



where Ht '= ^ X (0,T] with T > 0. We assume Dirichlet boundary conditions, 
namely O = I on <90. When f2 is the whole of K 3 we assume O — > I as \x\ — > oo, 
which is equivalent to assuming u* — >• at the boundary. Note that in certain 
situations appropriate radiation conditions have to be taken into account. 



3. Propagation of rotational waves 

3.1. Assumptions of the model. We are now discussing solutions to the fully 
nonlinear problem. As a first simplifying assumption, we consider a medium which 
can only experience rotation about one axis, the z-axis say. This means we can 
assume 

(0 1 0\ / cos ip — sin ip 0\ 

-10 0, O = sin 95 cos<yj 01, 
0/ V 1/ 
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where tp = (p(x,y, z,t) is an arbritrary function of the three spatial variables and 
time. This choice for u* means that u x = and u y = in Eq. (|2.3p and thus 
divtt = d x u z . To further simply the problem, we follow the well-known approach of 
classical elasticity [19] . Note the interesting twist waves in liquid crystals analysed 



Transversal rotational waves: We assume that the medium is homogeneous along 
the single axis of rotation, this means the z-axis. Thus, we restrict ourselves to 
ip = (f>t(x,y,t), which corresponds to the choice divtt = 0. 

Longitudinal rotational waves: By choosing ip to depend only on the z-coordinate 
and time, ip = <f>\(z,t), we effectively set curlv = 0. 

We will use this approach to find rotational waves propagating through the elastic 
medium. 

3.2. The Helmholtz equation. In general, since the matrix u* has three in- 
dependent components, the variation with respect to its components yields three 
coupled second order PDEs. Using the above ansatz for the transversal and lon- 
gitudinal rotational waves, two of the three equations are identically satisfied, and 
the remaining equations become 



respectively. Note that these equations are linear, however, we did not linearise. It 
is our ansatz that makes the nonlinear terms disappear. 

Therefore, in our elastic medium we find two types of rotational waves. Due to 
the nonlinear nature of the equations (f> = <fi t + (j>\ is not a solution of the Euler- 
Lagrange equations. As expected, the superposition principle does not hold. 

The transversal rotational waves travel in the x — y planes and are homogenous 
along the z-axis while the longitudinal rotational waves travel along the z-axis only. 
The corresponding wave velocities are given by 



It is a well known fact in classical elasticity that the wave velocity of the longitudinal 
waves is always greater than the transversal wave velocity. This does not hold in 
our rotational elasticity model. For an ordinary material we found c\ > 3c 2 + c 3 . 
Inserting this into (|3.2p the ratio of the velocities becomes 



which is in agreement with classical elasticity. On the other hand, for an auxentic 
material 2c 2 < c\ < 3c2 + C3, and thus 



in 0. 



(02 + c 3 )(d xx + d zz )(j)t - pd t t<t>t = transversal 
2 

-(ex + 2c 3 )<9 zz 0! - pdu<t>\ = longitudenal 




with their ratio given by 



(3.2) 
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These considerations show that our rotational elasticity model shares many features 
with well know classical elasticity, however, it also contains many new interesting 
features. 

3.3. Visualising transversal rotational waves assuming planar symmetry. 

If we are now separating variables and assume 

(3.3) 4> t (x, y, t) = cos(u>t)v(x, y) 

we obtain 



(3.4) (A + /c 2 )w = 0, k 

V c 2 + c 3 

the Helmholtz equation and the dispersion relation for the medium. The Hclmholtz 
equation is a well understood equation within the theory of PDEs. However, we 
will further simplify our model. 

Note that, as expected, the function <j> = const, solves these equations and hence 
it is natural to regard this solution as the ground state where all material points of 
the medium are aligned. 

In addition to assuming rotations about the z-axis only and homogeneity along 
the axis of rotation, we will now also specify a point in the plane orthogonal to the 
axis of rotation, the x — y plane. This means, we are assuming rotational symmetry 
in the x — y plane, we write <j)(x, y, z) = e lujt v(r) with r 2 = x 2 + y 2 . The nonlinear 
problem now reduces to a single ordinary differential equation which is given by 

(3.5) v"(r) + -v'(r) + P " v(r)=0, 

r c 2 +c 3 

and is solved by a linear combination of Bessel functions of first and second kind. 
We require the solution to regular at r = 0, therefore choosing one constant of 
integration to eliminate the Bessel function of second kind. We are left with the 
solution 



(3.6) v(r) =vqJo[u 

V V c 2 + c 3 

where vq — v(Q) and Jo denotes the Bessel function of first kind. 

This solution of rotational elasticity can now be visualised in the following man- 
ner. Let us attach to every material point in the medium an arrow to indicate its 
orientation. This is the practical realisation of attaching an orthonormal basis to 
every point in space. The choice cf> = corresponds to the situation where all ar- 
rows are parallel to the x-axis, which means that the medium under consideration 
is assumed to be aligned. We choose U!y/p/ {ca + C3) = 1 and uq = it, where the 
latter means that we rotate the central material point by a half turn. Thus we 
investigate the effect of the resulting rotational wave on the rest of the medium. 
The medium in its ground state <fi = in shown in the left panel of Fig. Q] while the 
distorted medium is shown in the right panel. 

4. Conclusions 

A new point of view was introduced to model rotational elasticity in a nonlinear 
setting using orthogonal matrices as the unknown variables. We showed how this 
theory can be related to classical elasticity and identified parameter ranges where 
the rotational medium would correspond to an auxetic material. After identifying 
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Figure 1. Rotational elasticity visualised. The left panel shows 
the undistorted medium where all material points are aligned par- 
allelly to the x-axis. The right panel depicts the solution of the 
transversal rotational wave equation when the central material 
point is rotated by a half turn. 



the most general energy functional of this model, two types of plane wave solutions 
were constructed analytically, they are solutions of the nonlinear Euler-Lagrange 
equations. These waves correspond to traversal rotational waves and to longitudinal 
rotational waves. It is interesting that in our rotational setting the transversal wave 
velocity can be greater than the longitudinal wave velocity. From a PDE point of 
view this model is very interesting since concepts like coercivity and convexity are 
difficult to establish due to the variables being matrices. It is therefore natural to 
consider next the question of existence and uniqueness of the Dirichlct boundary 
value problem. 
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